hist(airquality$Wind,xlab='Wind')

boxplot(airquality$Wind,xlab='Wind',ylab='Speed(mph)')
boxplot(Wind~Month,airquality,xlab='Month',ylab='Speed(mph)')

plot(airquality$Wind,airquality$Temp)
with(airquality,plot(Wind,Temp))
title(main='Wind and Temp in NYC')

with(airquality,plot(Wind,Temp,main='Wind and Temp in NYC',type='n'))

with(subset(airquality,Month==9),points(Wind,Temp,col='red'))

with(subset(airquality,Month%in%c(6,7,8)),points(Wind,Temp,col='blue'))
fit<-lm(Temp~Wind,airquality)
abline(fit,lwd=2)
legend('topright',pch=1,col=c('red','blue','black'),legend=c('Sep','May','Other'))

par(mfrow=c(1,2))
hist(airquality$Temp)
hist(airquality$Wind)
par(mfrow=c(1,1))
boxplot(airquality$Temp)

par(mfcol=c(2,1))
hist(airquality$Temp)
hist(airquality$Wind)

